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. Abstract. Non-relativistic three-dimensional magnetohydrodynamic simulations of Poynting-flux- 

^ ■ dominated (PFD) jets are presented. Our study focuses on the propagation of strongly magnetized 

Qj [ hypersonic but sub-Alfvenic flow (C^ << < V^) and the development of a current-driven 

<S ■ (CD) kink instability. This instability may be responsible for the "wiggled" structures seen in 

VLBI-scale AGN jets. In the present paper we investigate the nonlinear behavior of PFD jets in 
a variety of external ambient magnetized gas distributions, including those with density, pressure, 
and temperature gradients. Our numerical results show that PFD jets can develop kink distortions 
in the trans-Alfvenic flow case, even when the flow itself is still strongly magnetically dominated. 
In the nonlinear development of the instability, a non-axisymmetric mode grows on time scales of 
order the Alfven crossing time (in the jet frame) and proceeds to disrupt the kinematic and magnetic 
] structure of the jet. Because of a large scale poloidal magnetic field in the ambient medium, the 

■ growth of surface modes (i.e., MHD Kelvin-Helmholtz instabilities) is suppressed. The CD kink 
j—{ ' mode (m — I) grows faster than the other higher order modes (m > 1), driven in large part by the 

. radial component of the Lorentz force. 

o ■ 

INTRODUCTION 

o ■ 

' Magnetohydrodynamic (MHD) mechanisms are commonly considered to be the likely 

ci ■ model for accelerating the wind/outflow from Young Stellar Objects (YSOs), X-ray 

j> I binaries (.XRBs), Active Galactic Nuclei (AGN), Microquasars, and Quasars (QSOs) 

I [see, e.g.. Ill and references therein]. The toroidal (azimuthal) component of the magnetic 

■ field, generated by the twisting of the magnetic field threading an accretion disk, plays 
an important role in magnetic acceleration [2]. High resolution observations of AGN 
show morphological structures, such as "wiggles (kinks)" or "bends", not only on kpc 
scales, but also pc scales and smaller ||3i]. Such helical distortions might be caused either 
by plasma instabilities or by precession of the jet ejection axis due to the gravitational 
interaction of binary Black Holes (BBHs) [4], BH/disk, or galaxies. 

The magnetically driven jet model may be supported by recent observations. The 
rotation-measure (RM) distribution for the parsec-scale 3C 273 jet has a systematic 
gradient across the jet, and the projected magnetic field vector is systematically tilted 
from the jet central axis [5]. Similar RM gradients also occur in several BL Lac objects 
i6]. These results indicate the presence of a toroidal component of the field inside jet. 

Based on both theoretical aspects and observational results, it is natural to investigate 
plasma instabilities that might occur in these observed morphological structures. In the 
MHD model the toroidal field component provides the coUimation of the jet by magnetic 
tension or "hoop" stress. However, it is well-known that such a cylindrical plasma 



column with a helical magnetic configuration is subject to MHD instabilities. These are 
usually divided into pressure-driven instabilities (PDIs) and current-driven instabilities 
(GDIs) [see, e.g., 7]. GDIs arise because the toroidal field component of a magnetically 
driven jet is equivalent to a poloidal (axial) electric current, and such "current-carrying" 
jets are susceptible to MHD instabilities. Kelvin-Helmholtz instabilities (KHIs), which 
are driven by velocity gradients, must be also taken into account [see, e.g., 8]. Because 
these above instabilities would occur together in real astrophysical situations, it might 
difficult to separate their individual effects. 

In the recent past, much less theoretical attention has been paid by the astrophysical 
community to PDIs [9, 10] and GDIs ful[i2i,i2J than to KHIs. KHIs have been consid- 
ered for the past few decades by many theoretical or numerical investigators [for reviews. 



see 



13|,114j, and references therein]. In general, sub-Alfvenic jets (Vj < Va) are stable 



against KHIs, while, super- Alfvenic, but trans-fast-magnetosonic jets (Va < ^ ~ ^fm) 
are KH unstable. For super-fast-magnetosonic jets, if they have a considerably large 
fast-magnetosonic Mach number Mfm(= ^/^fm)» then the jets gradually become sta- 
ble again. However, there are several effects that can stabilize jets beyond the Alfven 
surface: (1) toroidal (rotational) velocity HI 5(1 and toroidal magnetic field [3 [13], (2) 
considerable mass entrainment which increases rj (= pj/Pext.) [18, 19], and (3) external 
magnetic fields [18] or magnetized winds |20, 21]. Also, as found by [22], jets can main- 
tain their stability by keeping the average Aflven Mach number (Ma = Vj/Va) within the 
jet to order unity by concentrating poloidal magnetic flux near the jet axis, giving it a 
reasonably stiff "backbone". 

Appl, Lery, & Baty [23] investigated cold super-fast-magnetosonic force-free jets 
(VpM < Vj). In the linear phase, the fastest growing GDI (m = 1) is nearly independent 
of the radial profile of the pitch. Lery, Baty, & Appl [24] then performed numerical 
simulations based on their linear analysis. In the early non-linear phase, nothing in their 
numerical results indicate a possible disruption of jets due to GDIs (m = and 1). To 
our knowledge, only two numerical studies of pure GDIs for astrophysical jets have been 
performed. Todo et al. f25"] applied their results to YSO jets, and Nakamura, Uchida, & 
Hirose [26] applied them to a large scale (kpc scale) wiggled structure of AGN jets. 

In the present paper, we report 3-D non-relativistic MHD simulations of PFD jets and 
investigate GDIs in decreasing (density, pressure, and temperature gradients) magnetized 
atmospheres. PFD jets should be applicable to small scale AGN jets, such as those on 
parsec scales and smaller. We believe that the helical morphological features (wiggles 
and bends) indicate the presence of GDIs in highly magnetically dominated jets. 



NUMERICAL MODELS 

The numerical approach in the present paper assumes non-relativistic compressible ideal 
MHD, with gravity neglected. We simulate the non-linear system of time-dependent 
MHD equations in a 3-D Gartesian coordinate system (x,j,z). All physical quantities 
are normalized with a unit scale length Lq, typical density po, typical velocity Vq, and 
combinations thereof, po and Vao are the initial value of the density and Alfven velocity 
at the origin of the computational domain. Lq = 2R^q is the initial jet diameter at z = 0. 



The initial distribution of the external ambient inter-galactic medium (IGM) consists 
of a stratified and magnetized atmosphere with possible gradients in density, pressure, 
and temperature. We adopt an initial current- (and therefore force-) free magnetic con- 
figuration, 7(= V X 5) =0, where B is the magnetic field and J the corresponding cur- 
rent density. Under the assumptions of flux and mass conservation, the density will be 
proportional to the local magnetic field as p °^ where a is a free parameter. For 
a — 2, the Alfven speed is constant throughout the computational domain, while for 
a < 2, the Alfven speed decreases with distance from origin and for a > 2 it increases. 
The gas pressure P is assumed to be polytropic, P p^, where F = 5/3 is the poly- 
tropic index. The atmosphere is artificially bound to prevent it from expanding under 
its own pressure gradient by introducing a "pseudo-gravitational" potential designed to 
hold on to the atmosphere without significantly impeding the advancing jet li27ll . In 
the lower "boundary zone" (Zmin < z < 0) the boundary velocity is specified for all 
time as vj = v,j,{r,z)^ + v^{r,z)z, where r = {x^ This represents a continuous 

cylindrical MHD inflow, powered by non-linear torsional Aflven waves (TAWs) into the 
"evolved" region (z > 0) of the computational domain. The total computational domain 
is taken to be \x\ < x^ax, \y\ < ymux, and Zmin <z< Zmax, where 

'''•max 5 ymax — 16 (32i?jo), 

Zmin — —1 (2^jo), and Zmax — 20 (40i?jo). The numbers of grid points in the simulations 
reported here are Nx x Ny x = 261 x 261 x 729, where the grid points are distributed 
non-uniformly in the x, y, and z directions. High-resolution 3-D computations were per- 
formed on a FUJITSU VPP 5000/32R (9.6 GFLOP/s peak speed per processor) and 
required about 6 hours on the 32-processor machine. 



NUMERICAL RESULTS 

We concentrate on the solution in the region, —2.25 <x,y< 2.25, and 0.0 < z < 18.0 
and discuss the dynamical behavior of a typical case: the growth of GDIs in a PFD jet, 
with parameter a = I. By and p decrease gradually along z-axis and tend to ~ z^^ for 
z > 1 . We set the plasma beta (/3o = 2C^q/yVIq = 10"^) at the origin, and the ratio Va/Cs 
increases gradually along the z-axis. Throughout the time evolution, a quasi-stationary 
PFD flow (Fexb/Fioi, ~ 0.9, where Fexb Poynting flux and Ftot. total energy flux) is 
injected into the "evolved" region from the lower "boundary" zone. Fig. 1 shows the 
time evolution of the density and the velocity field in the x — z plane for the typical PFD 
jet. In the early stage (t = 24.0), a strongly magnetized helical jet powered by TAWs 
advances into the decreasing poloidally magnetized IGM, remaining in radial force-free 
equilibrium (7 x 5 ~ 0). This is precisely the transverse structure predicted analytically 
by l[l6ll . The front of the TAW (a fast-mode MHD wave) is decelerated due to gradually 
decreasing Va; the accumulation of B^ occurs behind this wave front. The toroidal 
magnetic pressure gradient force d{B^/S7t)/dz becomes large and works effectively 
at this front, strongly compressing the external medium along this propagating TAW 
front. At later stages the front becomes super-fast magnetosonic and a bow shock (MHD 
fast-mode shock) is formed. The accumulation of B(j,, due to the gradually decreasing 
Va, causes a concentration of axial current near central (z) axis, and the distribution of 
magnetic pitch (\B(j) / rB^l) become large. The Lorentz force breaks the quasi-equilibrium 
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FIGURE 1. Time evolution of a typical PFD jet. Gray-scale images of the density distribution (logarith- 
mic scale) are shown with the poloidal velocity field in the x — z plane att = 24.0 (top), t = 32.0 (middle), 
and t = 40.0 (bottom). 



balance in the radial direction, and the axisymmetric PFD jet is disrupted by a GDI kink 
mode (m = I) (t = 32.0). At late times (t = 40.0, see Fig. 1 bottom and Fig. 2 right) a 
large-scale wiggled structure appears, while the MHD jet is still sub-Alfvenic V^et < 
(except at the front itself; see Fig. 2 left). This PFD jet will continue to be subject to 
the GDI, even thought it is magnetically dominated throughout the computational time. 
Fig. 3 shows the closed circulating current system, consisting of a current that is co- 
moving with the PFD jet (close to central axis) and a return current that flows outside 
the magnetized wind (cocoon). The external magnetized wind reduces the velocity shear 
between the jet itself and external medium, thereby suppressing KHIs even if the MHD 
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FIGURE 2. Distribution of velocities (left) and magnetic pitch (right) parallel to the z-axis 



jets become super- Alfvenic (i.e., a kinetic energy flux-dominated [KFD]). 



^=40.0 1:0 




FIGURE 3. Gray-scale image of the axial current density distribution in x — z plane. 



DISCUSSION 

There now is observational evidence that, in at least some extragalactic systems, the 
AGN jets are slowly coUimated over a length scale of pc (>> BH -i- disk system) [e.g., 
M87,|2S|].Ifthe MHD process is indeed the acceleration mechanism, the Alfven surface 
will be at a distance of order the disk's outer radius i?disk- The observations can set 
an upper limit on this of a few tens of pc (the radius of the ionized disk lEoll ). From 
theoretical considerations, the optically-thick disk is expected to be at about i?disk ~ 20 
pc L3flJ. So, an understanding of sub- Alfvenic to trans-Aflvenic jet flow may be needed 
in order to understand jet dynamics and morphology on sub-pc to pc scales. 



CONCLUSIONS 

We have investigated the non-linear behavior of PFD jets powered by TAWs in the 
extended stratified atmospheres, and find the disruption of PFD jets due to a GDI 
kink mode, leading to the formation of wiggles. Even if jets become super-Alfvenic 
during their propagation, GDI disruptions can still occur. Growing modes of MHD 
KHIs for super-Alfvenic flows are not seen in our parametric survey, because of external 
force-free magnetized helical winds surrounding the well-coUimated jet, suppressing the 



growth of KHIs. GDIs might be a possible origin of the wiggled morphology of AGN 
jets are slowly accelerated and coUimated on pc scales. 
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